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Abstract 

In this paper we address the relationship between zero temperature Glauber dynamics and the 
diffusion-annihilation problem in the free fermion case. We show that the well-known duality transforma- 
tion between the two problems can be formulated as a similarity transformation if one uses appropriate 
(toroidal) boundary conditions. This allow us to establish and clarify the precise nature of the relation- 
ship between the two models. In this way we obtain a one-to-one correspondence between observables 
and initial states in the two problems. A random initial state in Glauber dynamics is related to a short 
range correlated state in the annihilation problem. In particular the long-time behaviour of the density 
in this state is seen to depend on the initial conditions. Hence, we show that the presence of correlations 
in the initial state determine the dependence of the long time behaviour of the density on the initial 
conditions, even if such correlations are short-ranged. We also apply a field-theoretical method to the 
calculation of multi-time correlation functions in this initial state. 
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1 Introduction 

It has been known for a long time that there is a relation between Glauber dynamics p| and the symmetric 
diffusion problem in the presence of annihilation and deposition of pairs of particles for a certain choice 
of the diffusion constant ||, which corresponds to the case in which this problem can be solved using free 
fermions ||, Q|. Using this relation, Family and Amar [^| have computed the time evolution of the particle 
density in the transformed state of the annihilation problem that corresponds to random initial conditions 
in Glauber dynamics at T = 0. They have shown that their result only agrees with the previously known 
results by Spouge |6| if one starts with zero initial magnetization in the Glauber problem. In all other 
cases the long time behaviour of the density shows a dependence on the initial conditions, a result which 
differs from the well known universal behaviour valid for random initial states in the annihilation problem 
[Q. This rises the question of the correspondence between initial states in the Glauber and annihilation 
problems and more generally the relation between observables in the two systems. In this paper we show 
that the duality transformation between the two systems is really a similarity transformation, if one uses 
sector dependent, toroidal boundary conditions in the Glauber model (see below). We show that a random 
initial state in Glauber dynamics is mapped through this similarity transformation to a state with nearest 
neighbour correlations. Such state is translationally invariant, which allow us to recover the result by Family 
and Amar and also to compute higher order correlators using a field-theoretic technique. The relation 
between these quantities in the transformed state and their counterparts in the random initial state with 
density 1/2 is explored and some calculations are explicitly done. Finally, we emphasize the role of the 
correlations in the initial state for the failure of the universality hypothesis. 

The structure of this paper is as follows: in section 2, we present Glauber dynamics in terms of a quantum 
spin chain and the definition of a Temperley-Lieb algebra in terms of the spin 1/2 operators. The associated 
Hecke algebra permits the construction of the similarity transformation to the reaction-diffusion model and 
also determines the boundary conditions of the system. The correspondence between observables in Glauber 
and annihilation dynamics is also addressed. In section 3, we discuss the relation between the initial states in 
the two problems. In particular, we study the mapping between the random initial state in Glauber dynamics 
and a short-range correlated state in annihilation dynamics. Using these results and those of section 2 we 
establish a correspondence between correlation functions in the two problems. In section 4, we show that 
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the short range correlated state is a translationally invariant state and we show how to calculate multiple 
time correlation functions in this state using a field-theoretic technique. The density is explicitly computed 
and shown to agree with Family and Amar's result. We study the relation between two-point correlation 
functions in these states and the random initial state with density 1/2 and recover the zero time correlations 
as a special case. Finally, in section 5, we present our conclusions. 

For simplicity we discuss only the case of zero temperature Glauber dynamics although our results are 
easily generalised to other problems with little modification, e.g. Glauber dynamics at finite temperature 
or the model of generalised dynamics considered by Peschel and Emery [§) , which maps by a duality trans- 
formation to a model of diffusing particles with pair annihilation and creation away from the free fermion 
line. 



2 The transformation law for the operators 

It is well known ]3|, that certain reaction-diffusion systems provide physical realisations of Hecke algebras. 
The time evolution of these processes is described by a master-equation which can be conveniently written 
using an operator formalism, in which one assigns to each configuration of the system a state vector in an 
Hilbert space Jl(| . The probability distribution is then represented by a state vector obeying a "Schrodinger" 
equation dt \ $ ) = —H\ *S> ) with | } = J2 n P(n,i)\ n) where P(n,t) is the probability of finding confi- 
guration n at time t and the set of different | n) is supposed to be orthonormal and complete. The operator 
H is a linear and in general non-hermitian operator encoding the rules of the stochastic process. For some 
systems of interest this operator can be written |3[ ^) as a sum of generators of Hecke algebras. Given that 
there exists a relation between reaction-diffusion systems and Glauber dynamics Q it is natural to ask if the 
evolution operator for this system can also be written in terms of Hecke algebra generators and if so, what 
conclusions can one draw from it. In order to show this we define the following operators on a lattice of L 
sites 

L 

# ± = Yl ( X - e ^-! ) 0- - ( e 2j + e 2j _ 2 - 1 ) ) 

eij-i = \ (1 + ^ ) 1<3<L 

e 2j = \(l + o]g].,) l<i<i-l (1) 
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and eo = e 2 L which will be given below. The operators <rj, ctJ are Pauli spin matrices at site j. 
The set ej (1 < j < 2L — 1) forms a Temperley-Lieb algebra uM, characterized by the relations 



e i = e o 



(2) 



ej ei = eiej if | j — i | > 2 . 

In order to construct e 2 L explicitly we define the set of operators gj (1 < j < 2L — 1) and the duality 
operator D by (l^] 




( 3 ) 

V • j=l J 

where Yllti 1 <?j is the ordered product of the t/j's and X will be either 1 or cr£ in this paper. The important 
points about these two operators are that both commute with gj for 1 < j < 2L — 2 and that (#2L-i -^) 2 = 
(X52L-i) 2 - When X = 1, we will call the corresponding duality operator D + and when X = <r£, we will 
call it £>_. As we will see below the choice of the operator D is directly related with the different types 
of toroidal boundary conditions that were referred in the last section. The operators gj, together with the 
commutation relations with X, form an affine Hecke algebra associated with the Temperley-Lieb algebra 
given above (see [|[ [l^] and references therein) and one finds D ej — e 3+1 D for 1 < j < 2L — 2. Since 
gj = (1 — i) ej — 1, one finds from (||) that gj = 1. Also X 2 = 1, and we conclude that D is unitary and 
hence invertible. We define e-ih as 

e 2 L = De 2L _ 1 D- 1 . (4) 

The set of operators ey (1 < j < 2L — 1) and e-ih satisfies the relations of a periodic Temperley-Lieb 
algebra [[[2) with 2L generators which is defined by (||) together with similar relations for e 2 L (e2L+i = e-i, 
etc). For the particular choices of X given above, it can be shown that Jll| 

e 2L = \ + Cal&l) if X = 1 

(5) 

e 2 L = £(1 - Ca z L &l) if X = a z L 
where C = J\j=i &j • ^ we substitute the definitions of the e/s in (|l|) H^ 1 can be seen to be the generator 
of the time evolution for Glauber dynamics. The + (— ) sign stands when X = 1 (X = a z L ). This explains 
the use of the notation H^. The choice X = 1 (X = a^), corresponds to the operator H + (H~) which, 



when applied to eigenstates of C with eigenvalue 1 (—1), generates the time evolution for Glauber dynamics 
at T = with periodic boundary conditions ]l4| j. When applied to eigenstates of C with eigenvalue —1 (1) 
H + (H~) generates the time evolution for Glauber dynamics with anti-periodic boundary conditions. Both 
dynamics are stochastic. We see that one can indeed write the evolution operator for Glauber dynamics in 
terms of Hecke algebra generators. We now define the following similarity transformation 

V± = RD ± (6) 

where R = exp ^i^ J2^=i 1S a global rotation of -| around the y axis. From the relations between the 
Temperley-Lieb generators e J+1 = D ej D^ 1 it can be easily shown that 

V ± of V± 1 = of of +1 1 < j < L - 1 
V± of of +1 V^ 1 = of +1 1 < j < L - 1 

(7) 

V±olV^ = ±Q L olof 
V±Cdib\V± 1 = ±3f 

where Q L = RC R^ 1 = l\f=i From (0) it also follows that V± CV± 1 = ±Q L . Applying V+ to H+ 
and V- to H~ gives 

L-l 

ff± = V ± H ± V^ 1 = y £-(i-df6f +1 )(l--(&f + vf +l )) 
+ 5(1 - 5f5f)(l-5(df + 

+ 5(1 t QLbroi){\- l -{b\ + bi)) . (8) 

The operators H^, restricted to the subspaces Ql = 1 for H + and Ql = —1 for H~ , are equivalent to 
the Hamiltonian of the diffusion- annihilation problem with rates of diffusion 1/2 and rate of annihilation 
1 and periodic boundary conditions, which can be solved in terms of free fermions ||. The other cases 
correspond to non-stochastic processes. Hence, we obtain a rigorous formulation of the well known duality 
transformation between the two models Q. These results are summarized in Table 1. Notice that, although 
we have not used it explicitly, the similarity transformation preserves the relations (|^) so one can also 
represent in terms of Hecke algebra generators. However, one has used here an Hermitian quotient of 
the algebra which is different from the one used in references (§[ ^] and which allows the Hamiltonian for the 
diffusion-annihilation to be written as a linear combination of Hecke algebra generators. 
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Glauber Dynamics 


Annihilation Dynamics 


Operators 


Boundary Conditions 


Operators 


Boundary Conditions 


H+,6 = 1 


Periodic 


H+, Q L = 1 


Periodic 


H~, C = 1 


Antiperiodic 


H~, Q L = -1 


Periodic 


H+, = -1 


Antiperiodic 


H + , Ql = -1 


Non physical 


i?-, (7 = -1 


Periodic 


H~, Ql = 1 


Non physical 



Table 1: Summary of the relations between Glauber and annihilation dynamics. On the left hand side we 
have the Hamiltonian operator which is equivalent to the Glauber-Ising Hamiltonian and the sector of the 
Hilbert space where that equivalence holds, indicated by the eigenvalue of C. On the right hand side we 
have the tranformed Hamiltonian and the sector in which it is equivalent to the diffusion-annihilation process 
(given by the eigenvalue of Ql)- 

Finally, let us consider the action of V± in a single <rj operator. Expressing the operators gi in terms of 
Pauli spin matrices in (^) one can, using the commutation relations for these operators show that 

V±a*V ± 1 = -cr\&l...&l . (9) 

Using (||) , one obtains the following transformation law for a pair of o"| af k < I 

V±^ k afV ± 1 - {V±alV^){V ± afV^) 

= ffl&Z. ..<%...&;) = al +1 ...af . (10) 

Note that taking I = k + 1 or k = L, I = 1 we recover the equalities ([?]) concerning <t|, 1 . Equation (|l^) will 
be useful below when we derive equalities concerning correlation functions. We now proceed to study the 
effect of V± in the states of the theory. 

3 The transformation law for the states 

The Glauber-Ising Hamiltonian at T= has two ground states, the ferromagnetic states with all spins up 
or down. Since H + is only equivalent to it in the subspace of the states with 0=1 we have to find a linear 
combination of these two states that belongs to this subspace (for simplicity we will specialize in V+). This 
state is simply | $ ) = 5 (|T • • • t ) + 14 • • • I ))• It is normalized in the sense that ( s \ ^ ) = 1, where ( s \ is 
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the sum of all configurations of spins with weight one. This expresses the fact that, for a stochastic process, 
the sum of the probabilities of the different configurations accessible to the system (e.g. the configurations 
of spins in the Glauber problem) has to add to one. Also if H is the stochastic Hamiltonian representing 
the dynamics then it follows from conservation of probability that (s \ is the left ground state of H, i.e. 
( s | H = 0. One can easily check that C \ ) = | ^ }• This equality implies that 

V+\9) = V+&\#) = {V+CV+ 1 )V+\ *} = Q L V+\*) (11) 

where we have used the transformation law for C found above. Thus, the transformed state V+ | ^ ) belongs 
to the eigenspace with Ql = 1. In the annihilation language this corresponds to the sector with an even 
number of particles B. Also, since H + \ ty) = 0, one obtains 

V+H + \V) = (V+ H + V+ 1 )V+ | *) = H + V + \ty) = (12) 

and hence V+ \ 9 ) is a ground state of the annihilation Hamiltonian (in the subspace Ql = 1, H + is 
equivalent to it). The only ground state belonging to the subspace with an even number of particles is the 
vacuum | 0), i.e. the state with no particles. Hence, we conclude that V+ | \P ) oc | 0). The proportionality 
constant can be shown, using (^), to be equal to -^=(— l) i ~ 1 e l T( i_1 ) and can be absorbed in the definition 
of V+ . Following the same reasoning that led to ( |l2"| ) , and in light that (s \ C — ( s \ , one can similarly show 
that ( s | V+ l H + = and that ( s \ belongs to the even sector. On the same grounds of uniqueness this 
shows that this state is equal (up to a normalisation constant which we absorb in the definition of V^ 1 ) to 
the left ground state of the annihilation Hamiltonian with an even number of particles, i.e. (s \ { even ; which 
is the sum of all configurations with an even number of particles with weight one. We are using the same 
notation ( s \ for the left ground states of the two Hamiltonians because both describe stochastic processes 
and this is the usual convention. 

We have also found that H~ is equivalent to the Glauber- Ising Hamiltonian in the subspace C = — 1. 
Since V- C VZ 1 = — Ql this subspace is mapped to the even sector of the annihilation problem. But in this 
sector H~ (@) is not equivalent to a stochastic Hamiltonian. If one starts with the subspace (7=1 then In- 
effectively maps this sector to the odd sector where H~ is equivalent to the annihilation Hamiltonian, but 
in this case H ~ is not equivalent to the Glauber-Ising Hamiltonian with periodic boundary conditions, but 
to the Glauber-Ising Hamiltonian with anti-periodic boundary conditions, which is also a stochastic process. 
Using the same argument as above one can show that the ground state of this Hamiltonian is mapped to the 
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ground state of the odd sector of the annihilation Hamiltonian. This state is just a uniform superposition of 
states with one particle at each site of the lattice If we apply VZ 1 to this state we will therefore obtain 
the ground state of the Glauber- Ising Hamiltonian with anti-periodic boundary conditions (the observation 
about normalisation factors also applies here). Such state is a uniform superposition of 2L — 2 states with 
two domain walls, one at the boundary and one at each site of the lattice, plus the two states with all spins 
up or down. These two states are, due to the anti-periodic boundary conditions, the image of the state 
with one particle at the boundary in the diffusion-annihilation model. The previous discussion shows that 
the transformations V+ and V— restricted to the C = 1 subspace, map a stochastic process (Glauber-Ising 
dynamics with periodic or anti-periodic boundary conditions) to a stochastic process (annihilation dynamics), 
something which we also pointed out in Table 1. In the other cases the mapping is from a stochastic problem 
to a non-stochastic problem. Nevertheless, these mappings can be considered from a purely formal point 
of view and examples of similarity transformations between stochastic and non-stochastic processes have 
already been studied in the literature |L4|, |l5|. |Tg| . 

Here we shall concentrate in the stochastic-stochastic mapping given by V+. We will consider the study 
of time dependent correlation functions in uncorrelated random initial states evolving in time according to 
Glauber dynamics. The state 



= n 



' '1+771 1 — 777- 

+ o? 

2 2 J 



*> (13) 



corresponds to the superposition of two random initial states with initial magnetization m and —to. From 
the point of view of the calculation of correlation functions of an even number of o| operators the two states 
are equivalent and | 4>) belongs to the (7 = 1 subspace. Under the application of V+, | $ ) transforms to 

i*> = til 1 -^ + H^w+il |o) (14) 

3=1 L J 

where we have used (^) and the fact that V+ \ ^ ) = | ). This is an initial state with short range correlations 
and its form will play a crucial role in the determination of the correlation functions in the late time regime. 
Under V+ the multiple time correlation functions of an even number of <rj spins at times t\ , ti , etc. , transform 

as 

( * i (*o (*o (**) ^ (*»)••• ^ m i $ ) 
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- \ s I CT ii+i • ■ ■ CT i2 e CT i 3 +i • ■ ■ a ji e 

x . . . e -« + a| 2N _ 1+1 . . . a| 2jv e-" +t ~ | $ ) 

= (s \<r? 1+1 (h) . ..a] 2 (h)a? 3+1 (h) . ..«rj» • ..^_ 1+1 N • ■ .a| 2N (^)|$) (15) 

where we have used equations (||) and (|To|) . We suppose that j\ < j'2, J3 < J4, etc. Also one concludes that 
for any odd number of <rj operators, one has 

(s\*l(t 1 )...&l N+l (t2N+l)\$) = 0, (16) 

since <tJ anticommutes with C . 



4 The calculation of correlation functions 



Using (p_5|) we have, in particular, that 



<*i|(i - ^-i(*)^(*))i^> = (s\\{i-m))\®) = • (17) 



The operator |(1 — ajLi crj) checks for the existence of a domain wall at site j in the Glauber problem, 
i.e. it checks if at site j the spins cease to point up and start to point down or vice-versa. The operator 
2 (1 — &j ) checks for the existence of a particle at site j, i.e. if spin j is down in the annihilation problem 
This indeed corresponds to the well known duality transformation [^|, |5|. Similar relations hold for higher 
order correlation functions. Now we will use the fact that the diffusion annihilation Hamiltonian can be 
completely solved in terms of free fermions by means of the Jordan- Wigner transformation pjj]. Using it, 
one is able to write the spin raising and lowering operators s+, sj (s^ = i (<tJ ± i (ij )) in terms of creation 
and annihilation operators aj, cij for spinless fermions. It turns out that the calculation of a multiple time 
correlation function of the fij operators in state | $ ) {fij ~ o,j in the fermion language) can be reduced 
to the calculation of objects like ( s \ b Pl . . . b P2k $ ) (k < N) jl8) where b p is the Fourier transform of the 
annihilation operator a^. The momentum labels pj are half-odd integers between — L/2 + 1 and L/2 || [|. 
Also, one is able to express the state (|l4[) in terms of fermion operators. The result is 

l*> = llf 1 ^ + ^<^+il 10) = exp(/3^(^ +1 - 1)| |0) 



i=i 



exp [ j3 ( qt qt +1 + atoj+a + a] +1 hj + a j+1 aj - l) J 0) (18) 



where m — e~ 2 ^ . The first equality follows from the fact that ( <jj <rj +1 ) 2 = 1 and the second just follows 
from the rules of the Jordan- Wigner transformation. In terms of the momentum space operators we can 
write | $ } as 



$) =cxp ( 2(3^2 [cos fe) 



(blb p + bl p b- p ) + sin ^ (b P b- p + P_ p bl) -1] |0> .(19) 



The simplest approach to use if one wants to calculate the density or any equal-time correlators is given in B . 
For the calculation of multiple time correlators, we will follow a different route. We will look for operators 
cp, c p that diagonalize the quadratic form appearing in the exponent of (|l9|). This can be accomplished by 
means of a Bogoliubov transformation |l9| 

c p = cos(^) b p - Bin (20 bl p 

(20) 

4 = cos(^) S| - sin(^) S_ p 
and one gets 



$> = exp ( 2/3^(4 



ctc p - 1)1 |0) . (21) 

One then expands the exponential using the anti-commutation relations for the c p , c p . Re-expressing 
these operators in terms of w pl b p and applying them to the vacuum, one finally obtains 

| I) = e-^nt> + *p6lpSt]|0) (22) 

p>0 

where 7 P = cosh(2/3) — sinh(2/3) cos(^| £ ) and 5 p = sinh(2/3) sin(^2). The form ( p2| ) expresses the 
translation invariance of the state | $ ). One defines the following object 



z [Vp,V- P } = (0 I exp I ^(cot(^)fe p S_ p + i] p b p + ij^ p b- p ) J |$) (23) 

\P>0 / 

where the quantities r] p , r\- p are Grassmann variables anticommuting among themselves and with the b p 's. 
Their presence is necessary to make the terms in the exponential commute with each other. It can be shown 
[^8| that Z[rj p , ?y_ p ] is the generating function for the quantities (s \ b Pl . . . b P2k | $ ) that is 

( s | b Pl . . . b P2k | * > = d Vpi ... Z[ Vp , 77_ p ] |„ p=0 . (24) 



Substituting (22) in (E3J) we find, after a few algebraic manipulations, involving the anticommutation 



relations between the b p , b p ITS, the following expression for Z[rj p , ry_ p ] 
Z[v P ,V-p] = exp sinh(2/3) e" 2/3 sin V-pVp^J ■ 



(25) 
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Since W — In Z is a quadratic function in the 77's it immediately follows that a Wick's decomposition 



holds for the quantities ( s | b Pl . . . b P2k | $ ) mM. In particular, from (24) and (|25|) one has 



(s\b p ,b p \ *) = sinh(2/3) e" 2/5 sin ( HlZ. J <5 P: _ P < . (26) 



L 

Given that the expression for the space-dependent average density is |4j] 

<»j(*)> = X £ e^'W-t'-'+v)* cot (^) ( s I V S_p I $ ) (27) 
p,p' 

where e p = 1 — cos p, one obtains, substituting ( p6f ) in ( p7|) and taking the thermodynamic limit L — » 00, 
the following expression for the density of particles pit) at time < 

p{t) = \{1- m 2 )e- 2t (I Q (2t) + h(2t)) (28) 

where Io(2t), I\(2t) are the modified Bessel functions of order zero and one, and where we have used the 
identity m = e -2 ' 3 . This is the well-known expression obtained by Family and Amar [j5| who have also 
considered a random initial state in Glauber dynamics. They have shown that while the initial distribution 
of spins is uncorrelated the distribution of domain walls, i.e. the distribution of particles in the annihilation 
problem, is correlated. This correlated structure is found in the transformed state | (|l4]). For m = 0, 
the expression (^8|) is identical to the one found by Spouge || for an uncorrelated random initial state with 
initial density 1/2 in the annihilation problem. Indeed, this identity is more general. If m = then this 
means that we have to take the limit (3 — > 00. If we take such a limit in equation ( p2|) then one obtains the 
exact expression for a random initial state with density 1/2, projected over the even sector Q]. Therefore, 
our calculations provide a rigorous framework for the well-known duality transformation. In the long time 
limit t — > 00 one finds from ( |2^ ) the leading behaviour p{t) w 2v^r? ' ^ ^ c P en ^ s 011 ^ nc initial conditions (i.e. 
magnetization) ||. The amplitude of p(t) differs from the universal result found for uncorrelated random 
initial states 0, [l^] . Thus it is seen that the presence of correlations in the initial state breaks the universality 
of the amplitudes of correlation functions. 

Our results, namely the Wick decomposition, also allow the explicit computation of higher order corre- 
lation functions. Apart from a factor 2 sinh(2/3) e _2/3 the result ( p6[ ) is identical to the results obtained for 
the random initial state with density 1/2 |is| and, as stated above, it reduces to it when m = 0. Therefore 
one has 

(s\b Pl ...b P2k \$) = (l-e-^) k (s\b Pl ...b P2k \l/2r en . (29) 
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As an example of the above result let us consider the two-point correlation function ( s | hj(t) n&(i')| 
We use (^9|) with k — 2 and k = 1 because this object can be written as a linear combination of the 
terms ( s \b Pl b P2 b P3 b Pi | $ ) and ( s \ b Pl b P2 | $ ) iQ, which contribute with different powers of (1 — e~ 4/3 ). 
Separating such powers, we obtain 

( a | h, (t) h k (t')\ $) = (!- e-^fis | h 3 (t) n k (t')\ 1/2 ) even 



i_ r n - e" 4 ' 3 ) - (1 - e- 4 ^) 2 ] x 



2L 2 



(30) 



{l]p P e 2 S i(pi_P2)(: '~ fc) ~ (e -fi +e P2)*+( e -pi- e -P2)*'sin (^2.) cot (^g 1 ) 
+ Dp p e ^ (pi ~ P2)(3 ~ fc) ~ (e -pi +e P2) t -( e Pi- £ P2) t ' s in ( 2 ^ Pl ) cot (^zr)} 
In the thermodynamic limit L — > oo, we have for t' = t 

(s\h 3 (t)h k (t)\$) = (l-e-^) 2 (s\hj(t)fi k (t)\l/2) even 

+ e -^(l-e-^) x jV 2 *(/ (2<) + h(2t))5 jtk 

+ -k)- 9(k-j)]e- 2t (Ij-k-im - J,-_ fc+1 (2t))| (31) 

where we have used the integral representation of the modified Bessel functions Ij(2t) and is the 
Heaviside step function. The term 0(j — k) — 9(k — j) is the Fourier transform of cot (^) pl| . In particular, 
when t = 0, this reduces to 

+ \ e" 4/3 (l " e" 4/3 ) (* A fc+i + ) (32) 

where we have used the fact that there are no correlations in the state | l/2) e " e ™ at t = 0. One clearly 
recognizes in the first two terms the contribution of the unconnected part of the correlation function. But 
the third term indeed confirms that even at t — there are short range correlations. This term is zero when 
we take j3 — > oo and we just obtain the trivial result for the | 1/2 y ven state. 

If we use equation ( |30| ) to calculate the density-density correlation function in the thermodynamic limit 
we see that the second term of (^0|) will vanish. This leaves us with a term with an amplitude proportional 
to (1 — m 2 ) 2 . The ratio of this correlation function with the square of the density (p8|) is independent of m in 



agreement with the general results known from the renormalisation group approach (see for example 20 ). 
So despite the fact that the amplitudes of the various correlation functions are non-universal as emphasized 
above, their ratios obey the universality hypothesis. 
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The results discussed above show that this approach allows not only to recover the know results but also 
provides a convenient way to compute higher order correlation functions that can of course be translated 
back to the Glauber-Ising language. 



5 Conclusions 

We investigated the relation between the Glauber-Ising model at zero temperature and the diffusion-annihi- 
lation model in the free fermion case. We obtained the following new results 

(i) The duality transformation between the two models can be formulated as a similarity transformation 
if one uses an Hamiltonian with sector-dependent toroidal boundary conditions. The transformation laws for 
the operators are explicitly given. We also obtain the transformation laws for the states. This permits a one- 
to-one correspondence between a state of Glauber-Ising dynamics and a state of the diffusion-annihilation 
problem. In particular an uncorrclatcd random initial state in Glauber dynamics transforms to a state with 
short range correlations. 

(ii) Using the free fermion solution of the diffusion-annihilation problem we have computed the time de- 
pendent behaviour of the density and equal time two point correlation function in this short-range correlated 
state. For the density, we recover the results of the literature. We show that, surprisingly, the presence of 
correlations extending over only one lattice site in the initial state leads to a long time behaviour of the 
density dependent on the initial condition. The field theoretic approach we have used can be applied to 
the study of higher order correlation functions. Its use depends on the form of the initial state which 
reflects the translation invariance of this state. An initial state with the Wick's decomposition property was 
also discussed by Balboni et al. j2l). They have considered a continuous system described by a boson field 
theory. The initial state which they have chosen is characterized by the fact that the connected correlation 
functions of the density operator of order higher than two vanish. They found that for pure annihilation the 



amplitudes are also non- universal. In the case of the initial state (14) the higher order connected correlation 



functions of the density operator are non-zero, which shows that this state has a different structure. 
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